An improved algorithm for quantum separability and entanglement detection 
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Determining whether a quantum state is separable or entangled is a problem of fundamental 
importance in quantum information science. It has recently been shown that this problem is NP- 
hard. There is a highly inefficient 'basic algorithm' for solving the quantum separability problem 
which follows from the definition of a separable state. By exploiting specific properties of the set of 
separable states, we introduce a new classical algorithm that solves the problem significantly faster 
than the 'basic algorithm', allowing a feasible separability test where none previously existed e.g. in 
3-by-3-dimensional systems. Our algorithm also provides a novel tool in the experimental detection 
of entanglement. 



Entangled quantum states are interesting both from 
theoretical and practical points of view. Theoretically, 
entanglement is connected to the confounding issue of 
nonlocality. Practically, entangled states arc useful in 
quantum cryptography and other quantum information 
processing tasks (see and references therein). A mixed 
quantum state is defined as separable if and only if it can 
be written as a convex combination of pure separable 
states (and defined as entangled, otherwise). Solving the 
quantum separability problem simply means determining 
whether a given quantum state is entangled or separable. 
The problem comes in two flavors - one theoretical, and 
the other experimental. In this paper, we describe an al- 
gorithm for solving the quantum separability problem in 
the theoretical setting. We also describe the algorithm's 
utility in the experimental setting. 

We begin by introducing some notation and precisely 
defining the quantum separability problem. In what fol- 
lows, we are considering a bipartite quantum system of 
dimension M x N . Let Hj\/jv denote the vector space 
of all Hermitian operators acting on C M ® <C N . Noting 
that Mm,n is isomorphic to M. M N , it is endowed with 
the Euclidean inner-product (X, Y) = tr(XY), which in- 
duces the corresponding norm ||A|| = y / tr(X 2 ) and dis- 
tance measure \ \X — Let T>m,n C Hm,jv denote the 
set of all density operators. The set of bipartite separable 
quantum states Sm,n C T>m,N is defined as the convex 
hull of the separable pure states {|a)(a| <g> |/3)(/3|}, where 
| a) (|/?)) is a normalized vector in C M (C N ). An arbi- 
trary density matrix in T>m.n is parameterized by n — 1 
real variables, where 

n = M 2 N 2 , (1) 

and an arbitrary separable pure state is parameterized 
by 

k = 2(M + iV) - 4 (2) 
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real variables. Of course, in defining the separability 
problem, we cannot allow infinite precision, so we need to 
introduce a precision parameter d > 0. We are now ready 
to define the (quantum) separability problem as follows: 
Quantum Separability Problem. Given a density 
matrix p £ T>m,n o,nd a precision 6, assert either 

SEPARABLE: there exists a separable state a 
such that \\p — a\ \ < S; or 

ENTANGLED: there exists an entangled state r 
such that \ \p — t\ \ < S. 

The separability problem has been shown to be NP- 
hard Q, thus any devised test for separability is likely 
to require a number of computing resources that scales 
exponentially with M and N. There exist efficient "one- 
sided" tests for separability, where the output of some 
polynomial-time computable function of the matrix for p 
can indicate that p is certainly entangled 0, 0, B IE or 
certainly separable @, EJ, EJ , but not both. 

The experimental flavor of the separability problem 
can be defined as follows: Given many physical copies 
of a completely unknown quantum state p, determine 
whether p is separable. One may solve this problem by 
performing full state tomography in order to construct 
the density matrix for p to some precision 5, and then 
solve the theoretical separability problem. If there is 
some partial knowledge of p, then there are more options, 
such as testing for a violation of a specific Bell inequality 
[III IT^ I or invoking entanglement witnesses 0, Q] . In 
the case where MN < 6, the positive partial trans pos e 
(PPT) test HOI 

can be implemented physically |16lll7| . 
though currently this approach is not experimentally vi- 
able. 

The 'basic algorithm' that follows from the definition 
of a separable state is simply a straightforward search for 
a convex combination of separable pure states that gives 
the required density matrix within precision 5. Since any 
separable density operator in Sm,n can be written as 
a convex combination of n separable pure states [l^, a 
lower bound for the worst-case run time, ti, of this search 
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is given by 

h(n,S) > ^ (J^]) xpoly(n,log(l/<5)), (3) 

where Q$ is the number of pure separable states to pre- 
cision S. The first binomial factor is a lower bound for 
the number of combinations of n distinct pure separable 
states to precision S; the second is a lower bound for the 
number of probability distributions over the n states to 
precision 5. The algorithm that we present here has a 
worst-case run time, ti, with an upper bound of 

t 2 (n,6) < n s xpoly(n,log(l/5)). (4) 

To compare these run times it is not necessary to compute 
the exact value of 0,5, but note that for realistic values 
of 5, fig will be of order 0(2/S k ) 3> n, so t 2 will be 
significantly less than t\. Later it will be explained that 
equation can be improved further in practice with the 
use of global optimization routines, making the algorithm 
of practical use in the case where M and N are small 
(and S is not too small). Even for M = N = 3, there was 
previously no known better algorithm for the separability 
problem than the basic one described above. Note also 
that for MN < 6, where the PPT test is necessary and 
sufficient, our algorithm still offers its novel advantage in 
the experimental setting (as explained later). 

Before describing our alg orithm for separability, we 
note the following fact [l5|: A state p is entangled if 
and only if there exists an entanglement witness [jjj that 
detects it. An entanglement witness is any traceless op- 
erator A G M.m,n for which there exists a state p G T>m,n 
such that 

tr(Ar) < tr{Ap) Vcr G Sm,n- (5) 

This definition is slightly different from that used in the 
literature, however it substantially simplifies the descrip- 
tion of the algorithm. Recalling that Mm,n is isomorphic 
to R n , the above definition implies that for entangled p 
there exists a hypcrplanc which separates p from the set 
of all separable states 5m, n- If one defines the function 

bA = max tr(Aa), (6) 

then the set {X £ Mm,n : tr(AX) — bA} is one such 
hypcrplane. We (non-uniquely) define a a to be any ele- 
ment of Sm,n such that tr (.Acta) = bA- It suffices only to 
consider entanglement witnesses A such that tr(A 2 ) = 1. 
that is, those which lie on the (n — 2)-dimensional sur- 
face of a O-centered unit-hypersphere in Hm,jv, where 
is the origin (null operator in Mm.n)- For our pur- 
poses, however, it will be useful to characterize all poten- 
tial entanglement witnesses by the corresponding (to — 1)- 
dimcnsional unit-hyperball, VV, defined as 

W = {A e M m ,n ■ tr(A) = 0, tr{A 2 ) < I}. (7) 



Entanglement Witness Problem. Given a density 
matrix p G X?m,jv o,ud a precision 5, either assert 

SEPARABLE: there exists a separable state a 

such that \\p — o~\\ < 5; or return 

A G W: an operator such that 

tr(Ao-) < tr(Ap) + S for all a G S m ,n- 

The witness problem is thus to decide that p is almost 
separable, or to find an approximate entanglement wit- 
ness for p. Note that any algorithm solving the witness 
problem also solves the separability problem. Our algo- 
rithm actually solves the witness problem. 

Our algorithm is an iterative one, which calls a com- 
putationally expensive subroutine at each iteration. It 
is convenient to treat this subroutine as a black box, or 
oracle, when describing the algorithm's main structure. 
Simply define the oracle, O, such that it takes an oper- 
ator A, and returns O(A) = a a- Whichever way O(A) 
is computed, it suffices that the maximization in © is 
done over the pure separable states (and that a a is a pure 
state). Thus, the motivation behind reducing the sepa- 
rability problem to the oracle O is quite simple: we are 
exploiting the fact that the separable pure states, which 
are the extreme points of 5m, jv, are parameterized by 
k variables rather than n. Thus, from a practical point 
of view, the complexity of computing O(A) scales much 
better than that of either a brute-force search through 
all entanglement witnesses, or the 'basic algorithm'. 

Before delving into the details, we give a high-level de- 
scription of our algorithm. The algorithm maintains a 
set K, C W of operators which arc potential entangle- 
ment witnesses for p. If (and only if) p is entangled, 
there exists a closed, convex subset of W, which we call 
W p , consisting of all entanglement witnesses that detect 
p. Throughout the algorithm, we have W p C K,. Initially, 
JC is set equal to W. In each iteration, the algorithm se- 
lects a test-witness, A G JC, and computes a a = 0(A). 
If tr(AcrA) < ti(Ap) + 5, then the algorithm returns A; 
else, /C is reduced and the next iteration begins. The 
algorithm keeps reducing the set K, until it either finds 
an (approximate) element of W p , or it decides that W p 
is empty, and therefore p is separable. To decide W p is 
empty means that tC is too small to contain W p . This 
requires having a lower bound on the size of W p . By ex- 
ploiting the role of d in the problem definitions, such a 
lower bound can be derived in terms of S and n. In what 
follows, we will ignore the role of 5, as it obfuscates the 
main idea of the algorithm. 

We now describe how to reduce the set fC, that is, 
to discard elements of JC that are not elements of W p . 
Suppose A is not in W p , but is sufficiently close to W p . 
Then, A, p, and a a can be used to define a half-space 
{X G M m ,n ■ tr(KX) > 0} that contains W p . Specifi- 
cally, we have the following: Let W be any operator in 
W p and suppose A £ W p . If tr(WA) > 0, then choosing 



We can now define the (entanglement) witness problem, 
a problem slightly harder than the separability problem: 
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(and then normalizing K) gives to(KA) = by construc- 
tion, and it is easy to verify that to(KW) > 0. The idea 
is that, at each iteration, the test- witness A is chosen so 
that it is (approximately) in the center of the current fC 
(relative to the Euclidean geometry). If the oracle O re- 
turns (7 a such that to{Ap) < to {Act a), then, as long as 
tr(WA) > for all W G W p , equation (jSJ) gives a cutting 
plane {X G Mm,n '■ to(KX) = 0} that slices through 
A and 0. This allows us to discard the half of K. con- 
sisting of operators X such that to{KX) < 0. Because 
}C is being approximately halved at each step, the algo- 
rithm quickly either finds an entanglement witness for p 
or concludes that p is separable. 

Our problem of determining whether the convex set 
W p is empty using cutting planes is well studied in the 
field of convex optimization. However, because of our 
special requirement that tr(WM.) > for all W G W p , 
none of the existing algorithms can be applied directly. 
Fortunately, though, the analytic-central-section algo- 
rithm due to Atkinson and Vaidya [2(j can be adapted 
for our purpose, giving an algorithm with the desired 
complexity. 

We now describe the algorithm. Let Imn be the 
maximally mixed state, which is properly contained in 
Sm,n P- llflj. It is easy to verify that W p must be con- 
tained in the half-space {X : to((p — Imn)X) > 0}. Let 
K\ = (p — Imn)/\\p — Jmjv||- Thus, straight away, fC is 
reduced to the half-ball W C\ {X : tr(A'i^) > 0}. The 
first test- witness to give to the oracle is A = p — Imn 
(which is along the center-line of the half-ball). If the 
oracle confirms that A detects p, then we are done. Oth- 
erwise, we use equation IjSJ to generate a cutting plane. 
By way of mathematical induction, assume that, at some 
later stage in the algorithm, JC has been reduced to 

K = Wp|nti{* : to{KiX) > 0}, (9) 

by the generation of h cutting planes {X : to(KiX) = 0}, 
as described above. Recall that we want to choose a test- 
witness that is approximately in the center of JC. An 
easily computable candidate is the analytic center, O, of 
JC |2l| . which is defined as the unique minimizer of the 
real convex function 

h 

F(X) ee -^log(tr(^A))-log(l-||A|| 2 ), (10) 

i=l 

defined for X G fC. The relation X7F(C) = gives 
q _ i- 1 \C\ | ^{KiC) ' w hi cn > by the inductive hy- 

pothesis, implies that tr(WC) > for all W G W p . Thus, 
A = C is a suitable test-witness to give to the oracle and 
to use in equation JSJ|. Full details of a robust algorithm 
are too numerous to include here but can be derived with 
the help of [20l l2ll 12^ . The important point is that the 
separability of a given density matrix can be decided with 
only n x polylog(n, 1/5) calls to the oracle. 

Now consider the complexity of computing 0(A), 
which so far has been black-boxed. The most naive way 



to carry out this computation is to one-by-one calculate 
to (Act) for each of the fig pure separable states a (to 
precision <5) and return the a that produced the largest 
value of to(Aa). Even with this naive way of computing 
O(A), the total run time of our algorithm is significantly 
shorter than that of the 'basic algorithm' for quantum 
separability (compare and Q). However, for any 
given orthogonal Hermitian basis of Mm,n, the closed, 
general form of the function tr(^4er) can be written down 
in terms of the k real parameters of the separable pure 
states. Armed with the closed form of the function to 
be maximized, various well-studied global maximization 
techniques are at one's disposal, for example, Lipschitz 
optimization or interval analysis |24| ■ Call the func- 
tion to be maximized / and denote its global maximum 
by /*. As the global optimization algorithm proceeds, it 
gives progressively better lower and upper bounds on /*. 
Call these bounds / and /, respectively. A key advan- 
tage of our algorithm is that, during any computation of 
O(A), the search for /* may be halted early when either 
(?) to(Ap) < /, in which case equation (JSJ) can be invoked 
to generate a new cutting plane, or (ii) f < to(Ap), in 
which case the algorithm has found an entanglement wit- 
ness for p. Thus, the algorithm's run time may be signif- 
icantly shorter than the worst-case analysis predicts. 

Finally, we discuss how the algorithm may be used 
when only partial information about the state p is avail- 
able. This is of particular use in an experimental setting. 
Let B be an orthonormal, Hermitian basis for Hjv/,jv- The 
state p can be written p = X)j=i PXiXi, where px { G M. 
Each coefficient px { is simply the expected value of Xi, 
which equals to(Xip). The expected values of all elements 
of B constitute complete information about p. Suppose 
we have only measured j < n expected values. The al- 
gorithm can be applied in this reduced, j-dimcnsional 
space. If the algorithm finds a hyperplane separating 
p from Sm,Ni then p is entangled; otherwise p may be 
entangled or separable, as the j expected values are con- 
sistent with a separable state. As expected values are 
being gathered through experimental observation, they 
may be input to the algorithm. If the basis B is sepa- 
rable, then the entire procedure can be done when the 
subsystems are spatially separated with local operations 
and classical communication. The idea of searching for 
an entanglement witness in the span of operators whose 
expected values are known was discovered independently 
and applied, in a special case, to quantum cryptographic 
protocols in pEj . 

We have given a classical algorithm for the quantum 
separability problem which takes as input the density 
matrix for a quantum state p and either decides that p 
is separable, or returns an entanglement witness that de- 
tects p. Our algorithm is the best-known algorithm for 
the general separability problem; it gets its advantage 
over the 'basic algorithm' from reducing the problem to 
an optimization over the pure separable states. If prop- 
erly implemented, the algorithm should give a feasible 
test for separability in low dimensions (e.g. MN < 10, 
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with current technology) . The general technique depends 
only on the convexity of the set of separable states, and 
thus can, in principle, be applied to test for multi-partite 
entanglement. The algorithm also gives experimentalists 
a tool for potentially determining if an unknown state 
is entangled by measuring only a subset of the expected 
values which completely describe the state. This method 
effectively trades quantum resources (additional copies 



of p) for classical resources (a computer able to calculate 
O). 

We would like to thank Carolina Moura Alves, Coralia 
Cartis, and Tom Stace for useful discussions. We ac- 
knowledge support from the EC under project RESQ 
(IST-2001-37559). LMI, BCT, and DC also acknowledge 
support from, respectively, CESG and NSERC; CMI; and 
NSERC and the University of Waterloo. 



[1] M. A. Nielsen and I. L. Chuang, Quantum Computation 
and Quantum Information (Cambridge University Press, 
Cambridge, 2000). 

[2] L. Gurvits, in Proceedings of the thirty-fifth ACM sym- 
posium on Theory of computing (ACM Press, New York, 
2003), pp. 10-19. 

[3] A. Peres, Phys. Rev. Lett. 77, 1413 (1996). 

[4] M. Horodecki and P. Horodecki, Phys. Rev. A 59, 4206 
(1999). 

[5] M. A. Nielsen and J. Kempe, Phys. Rev. Lett. 86, 5184 

(2001) . 

[6] A. C. Doherty, P. A. Parrilo, and F. M. Spedalieri, Phys. 

Rev. Lett. 88, 187904 (2002). 
[7] A. C. Doherty, P. A. Parrilo, and F. M. Spedalieri, A 

complete family of separability criteria (2003), quant- 

ph/0308032. 

[8] S. L. Braunstein, C. M. Caves, R. Jozsa, N. Linden, 
S. Popescu, and R. Schack, Phys. Rev. Lett. 83, 1054 
(1999). 

[9] K. Zyczkowski, P. Horodecki, A. Sanpera, and M. Lewen- 
stein, Phys.Rev. A 58, 883 (1998). 
[10] L. Gurvits and H. Barnum, Phys. Rev. A 66, 062311 

(2002) . 

[11] J. S. Bell, Physics 1, 195 (1964). 

[12] J. F. Clauser, M. A. Home, A. Shimony, and R. A. Holt, 

Phys. Rev. Lett. 23, 881 (1969). 
[13] O. Giihne, P. Hyllus, D. Brufi, A. Ekert, M. Lewen- 

stein, C. Macchiavello, and A. Sanpera, Phys. Rev. A 

66, 062305 (2002). 
[14] M. Barbieri, F. D. Martini, G. D. Nepi, P. Mataloni, 



G. M. D'Ariano, and C. Macchiavello, Experimental de- 
tection of entanglement with polarized photons (2003), 
quant-ph/0307003. 

[15] M. Horodecki, P. Horodecki, and R. Horodecki, Phys. 
Lett. A 223, 1 (1996). 

[16] P. Horodecki and A. Ekert, Direct detection of quantum 
entanglement (2001), quant-ph/0111064. 

[17] H. Carteret, Noiseless circuits for the Peres criterion 
(2003), quant-ph/0309216. 

[18] P. Horodecki, Phys. Lett. A 232, 333 (1997). 

[19] B. M. Terhal, Phys. Lett. A 271, 319 (2000). 

[20] D. S. Atkinson and P. M. Vaidya, Mathematical Pro- 
gramming 69, 1 (1995). 

[21] Y. Nesterov and A. Nemirovskii, Interior- Point Poly- 
nomial Algorithms in Convex Programming (SIAM, 
Philadelphia, 1994). 

[22] J. Renegar, A Mathematical View of Interior-Point 
Methods in Convex Optimization (MPS-SIAM, Philadel- 
phia, 2001). 

[23] R. Horst and P. Pardalos, eds., Handbook of Global 
Optimization (Kluwer Academic Publishers, Dordrecht, 
1995). 

[24] E. Hansen and G. Walster, Global Optimization Using 
Interval Analysis (Marcel Dekker Incorporated, Boston, 
2004), ISBN 0824740599. 

[25] M. Curty, M. Lewenstein, and N. Liitkenhaus, Entangle- 
ment as precondition for secure quantum key distribution 
(2003), quant-ph/0307151. 



